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This article presents the derivation of the stress-energy tensor of a free scalar field with a general 
non-linear dispersion relation in curved spacetime. This dispersion relation is used as a phenomel- 
ogical description of the short distance structure of spacetime following the conventional approach 
' of trans-Planckian modes in black hole physics and in cosmology. This stress-energy tensor is then 

used to discuss both the equation of state of trans-Planckian modes in cosmology and the magni- 
£Nj ' tude of their backreaction during inflation. It is shown that gravitational waves of trans-Planckian 

momenta but subhorizon frequencies cannot account for the form of cosmic vacuum energy density 
observed at present, contrary to a recent claim. The backreaction effects during inflation are con- 
firmed to be important and generic for those dispersion relations that are liable to induce changes in 
the power spectrum of metric fluctuations. Finally, it is shown that in pure de Sitter inflation there 
is no modification of the power spectrum except for a possible magnification of its overall amplitude 
' independently of the dispersion relation. 

> ■ 

qq . PACS numbers: 98.80.Cq 

<N 
On 

o ' !• introduction 

The inflationary paradigm provides an appealing framework to describe the very early phase of the evolution of 
. the Universe, notably because it produces in a natural way the seeds necessary to the formation of the cosmological 
large scale structures {!]. These initial density fluctuations are generated during inflation with an almost scale 
invariant power spectrum through the parametric amplification of the quantum fluctuations of the inflaton scalar 
field, with the accelerated Friedman- Lemaitre scale factor acting as a classical pump field. This can also be seen as 
particle production due to the breakdown of adiabaticity in the evolution of the quantum modes of the field as the 
wavelengths are stretched beyond the horizon. 

In spite of their successes, most models of inflation are subject to the trans-Planckian problem, namely the phase of 
accelerated expansion lasts sufficiently long that cosmological length scales today corresponded to scales much smaller 
than the Planck scale at the beginning of inflation ^jfy. Depending on one's point of view, this can be seen either 
as a problem, i.e. the celebrated predictions depend on unknown trans-Planckian physics, or as a blessing, since 
inflation then opens a window on physics beyond the Planck scale. Several studies have recently tackled the issue 
of the robustness of inflation to changes in super-Planck physics by adopting a phenomenological approach initially 
developed in the context of black hole physics, where a very similar problem arises |^,^), and in close relation to 
analogous problems in condensed matter physics [fjj. It consists in modifying the standard dispersion relation of a 
free scalar field for wavelengths smaller than the Planck length and in calculating the resulting power spectrum of 
inflation produced metric fluctuations ^,^,^,|). This approach is further motivated by the fact that the evolution 
of the scalar and tensor modes of metric fluctuations can be adequately described by free scalar fields propagating 
in a background spacetime ||. It should be noted that this is only a phenomenological approach to the problem of 
trans-Planckian physics, and that it is not unique. For instance, modifications of the canonical quantum operator 
commutation relations inspired from string theories [ p)[ have also been shown to affect the power spectrum of metric 
fluctuations 

Interestingly, concrete examples of dispersion relations leading to modifications of the power spectrum with respect 
to the standard (linear dispersion relation) predictions have been exhibited |^,^] but only in the case where the 
dispersion relation becomes complex. Even though complex dispersion relations are ordinary in classical physics or 
in quantum mechanics, they represent a problematic situation in the context of quantum field theory. It has also 
been shown that dispersion relations such that the evolution of the quantum mode is adiabatic all throughout the 
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inflationary phase up to horizon crossing cannot lead to significant modifications of the power spectrum || . With 
respect to the above example of dispersion relations, indeed adiabaticity is broken at the point where the physical 
frequency vanishes and the dispersion relation becomes complex. However it has been argued that such dispersion 
relations which break adiabaticity would lead to a possibly severe backreaction problem, i.e. the energy density 
contained in the modes would become greater than the background energy density ]l5| , [l6[ | . Therefore modifications to 
the power spectrum may arise if adiabaticity is broken, but at the expense of the creation of a possibly large amount 
of energy density. Unfortunately this amount could not be quantified rigorously as the authors did not have at their 
disposal the stress-energy tensor of a theory with modified dispersion relation. 

In this context, our main objective is to present a rigorous derivation of the stress-energy tensor of a free scalar 
field with non-linear dispersion relation. Such a dispersion relation breaks local Lorentz invariance as it implies 
the existence of a preferred reference frame, however it is crucial to maintain general covariance in order to achieve 
consistency with the Einstein equations (and notably the conservation of the stress-energy tensor). This problem 
has been examined in a number of studies, where an effective general covariant Lagrangian with explicit Lorentz 
invariance breaking has been constructed with the introduction of a dynamical unit timelike vector field whose role 
is to define the preferred rest frame . An effective Lagrangian describing a free scalar theory with a quartic 

dispersion relation could then be constructed along these lines, and the stress-energy tensor has been calculated for 
the particular case of the Corley-Jacobson dispersion relation Jl^Jl9t] , In the present paper we present a non-trivial 
extension of this latter study to the general case in which the squared frequency is a general analytic function of the 
squared momentum. Note that in Friedman-Lemaitre-Robertson- Walker (FLRW) cosmology, a preferred rest frame 
exists and coincides with the homogeneous isotropic spatial sections. Finally one should point out that the method 
given below cannot be used to derive the stress-energy tensor for theories with unmodified dispersion relations but 
modified canonical commutation relations as considered in Refs. [p^!3[. 

We then use this stress-energy tensor to address two points raised recently in the literature. We first discuss the 
claim p0| according to which a bath of gravitons of super-Planck momenta but frequencies much smaller than the 
Hubble expansion rate (hence a particular non-linear dispersion relation), can explain the form of vacuum energy seen 
in the Universe today. In order to do so, we calculate the energy density, the pressure and the equation of state of 
these quanta; we show that these gravitons possess neither the correct energy density nor the correct equation of state. 
As a second application we discuss the issue of backreaction of trans-Planckian modes in inflationary cosmology, using 
the expression for the energy density contained in trans-Planckian modes. 

This article is organized as follows. In Section || we formulate a covariant Lagrangian (§ II A) including extra terms 



to provide a modified dispersion relations an d calc ulate the corresponding stress-energy tensor (§ [II B| ), Then, we 
specify our result to the case of cosmology (§ |ll C|) . Section HI discu sses the two main applications of o ur res ults, 



n amely the equation of state of the trans-Planckian modes (§ III A) and the backreaction problem (§ IIIB and 



§ III C). Our results are summarized in Section IV. For the sake of clarity, we gathered notations and derivations of 
various identities used in the calculation of the stress-energy tensor in Appendix [A|; Appendix [B] presents the detailed 
derivation of the stress-energy tensor in the simpler case of the Corley-Jacobson dispersion relation. We use natural 
units in which h = fee = c = 1, and the metric g^ v carries signature (—,+,+,+). 

II. COVARIANT LAGRANGIAN AND STRESS— ENERGY TENSOR 

In this section, we first introduce a general covariant formulation of a Lagrangian describing a free scalar field 
with modified dispersion relation following the procedure described by Jacobson and Mattingly [ ^8|Jl9| , and derive 
the corresponding energy-momentum tensor. We try to remain as general as possible with r espe ct to the background 



spacetime M. and its metric, and defer the detailed study of a FLRW spacetime to Section UC. 



A. Definitions and covariant Lagrangian 

The action for a free scalar field with modified dispersion relation takes the form |p^| , ^9| 

5 t = |d 4 ^,+c+a (i) 

where C. is the standard Lagrangian of a minimally coupled free scalar field 



C^~-g^d^d^. (2) 
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Two contributions have been added to this Lagrangian in order to introduce the modified dispersion relation. As 
already mentioned in the introduction, a modified dispersion relation breaks local Lorentz invariance. In such a 
situation, a covariant formulation of the corresponding theory can be carried out by introducing a unit timelike vector 
field m m defining a preferred rest frame JlJ,|l^l . In Eq. (|l|) , the first term £ cor is responsible for the non-linear part of 
the dispersion relation and C u describes the dynamics of the vector field u^. These two corrective Lagrangians take 
the form 

£„, =-J2 b w ( V2n ^ ( p2 ^) ' £ " = -\(g^u^ + 1) - dxF^F^. (3) 

n,p<n 

The tensor F^ v is defined by F^ v = V M u v — V„ w M , where is the covariant derivative associated with the metric 
<7 M „; A is a Lagrange multiplier and the coefficients b np and d\ are arbitrary. The derivative operators T> 2n are defined 
further below. The overall Lagrangian maintains general covariance, which will notably ensure the conservation of 
the stress-energy tensor. The value of the Lagrange multiplier A can be obtained by the extremization of the action 
with respect to the vector field u M . 

Let us define more precisely the quantities appearing in the two extra Lagrangians in Eq. (|). We first assume 
that the spacetime M. is globally hyperbolic so that it can be foliated as M. — £ x R where £ are three dimensional 
spacelike hypersurfaces of constant q, where q is a scalar. It follows that the unit timelike vector field normal to 
these hypersurfaces is = — (d u q) / {~g a gd a qd^ q) 1 ^ 2 which indeed satisfies u M it M = —1, with respect to which we 
will define "time" and "space" components pf|]22] |. Frobenius theorem |2^] guarantees that the vector field u M is 
rotation- free and thus that the field strength tensor can be written as F^ — a^u v — a„ u M , a M being the acceleration 
defined in Appendix |A[ As a consequence, if the vector field is geodesic then F^ = 0. This will be for instance 
the case in the FLRW case but docs not need to be true in general. The projector on the hypersurfaces £ defined by 

-L/«/= 9^v + (4) 

coincides with the spatial metric as defined by an observer comoving with since the line element can be rewritten 
as 



ds 2 = g^dx^dx" = ~{u t ,dx l ") 2 + _L^„ dx^dx". (5) 

The covariant derivative T> a associated with the induced 3-metric of any tensor field Tp\'.'.'!p n , appearing in Eq. (|^) 
defining £ cor , is defined as 

^aVl— Ml ' ■ ' - L M„- L o'i ■ ■ ■ J -(J P - L a V /3 J Ai...A„ ' W 

By construction, it is the covariant derivative associated with Jl^ v and is orthogonal to u^, i.e. 

Various identities satisfied by u^, Jl^ v and D„, which are used repeatedly in the rest of this section are given in 
Appendix |a|. We finally define the operator V 2n appearing in the Lagrangian Eq. (||) as 

V 2n ee V^V^ 1 . . . V^W" . (8) 

For the particular case of a scalar field, this double derivative can be written as 

V^4> =±» V Q ±g =± af3 V a Vp<j> + u a V a <j>VpvP. (9) 

Throughout the paper, it will be understood that a derivative operator applies directly and only on the first term 
appearing on its right; derivatives of an ensemble of terms will be indicated using brackets. We will also use the 
shorthand notation of an overdot for the time derivative defined by an observer comoving with u' 1 , i.e. T^\'.\ M ^ n = 
u a V a T^l'.'.'.^ p n . In particular, <f> = u a W a 4>. The above double derivative operato r ca n be shown to coincide with the 
three-dimensional Laplacian as defined by an observer comoving with it M [see Eq. ( |A5| )]. The corrective Lagrangian C cot 
thus contains only "spatial" derivatives. For example, in Minkowski spacetime, this would yield plane wave solutions 
to the field equations with a dispersion relation for the pulsation lo 2 as a series in powers of squared momentum k 2 , 
where the first term in k 2 results from the free Lagrangian and the higher order terms come directly from the higher 
order Laplacians in £ cor . 
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B. Stress— energy tensor 



The stress-energy tensor is obtained by varying the action ([!]) with respect to the metric 



1 



- / d^xV^T.Jg^ = - / d^x^T^Sg^ 



(10) 



In the following, we derive this stress-energy tensor and extract the pressure and energy density under some hypotheses 
on the derivative D^. In Appendix [b], we give a detailed calculation of the stress-energy tensor in the simplest non- 
trivial case in which only b\\ is non- vanishing, as considered by Jacobson and Mattingly |l8| , |l9| . The covariant 
derivative obviously satisfies \T> 2 , g^,v\ = and we further assume that it also satisfies 



[V 2 ,u»] = [£> 2 ,V M ] =0. 



(11) 



Although these requirements seem restrictive and may not apply in general, they are fulfilled for the relevant cosmo- 
logical case of the FLRW spacetime we are interested in. We emphasize that these requirements are not necessary to 
the derivation of the stress-energy tensor in the case discussed in Appendix in which only the first coefficient b n 
is non- vanishing. The main point of the calculation is the variation of the corrective term S col: to the action, which 
gives 



E 5 « 



i=0 



v 2p 4> + v 2, \ 



■p-1 

.i=0 



-gd x. 



(12) 



As a result of the commutation relations, the second integral can be rewritten through multiple integrations by parts 
as 



£>2(n+p-l) ( 



5V A 



gd 4 x = / E{(j))5V 2 



-gd x, 



(13) 



where we used the short-hand notation E{<f>) = b np (n + p)V 2n+2p 2 <fi. Using expressions (B3-B4), we finally obtain 
the stress-energy tensor 



T =8 



+2 



fiu + 2Xu fl u L , 



n.p 



(14) 



-2V {ll EV v) </) - 2Eu [p y u) 4> - 24>u { ^J u) E+ V a E\7 a <j) + g^E<p + T$ , 

where T^) = —^diF^\F w x + g^diF" 13 F a p- The expansion 6 =±^ v V M u„ is defined in Appendix |A| together with 
the acceleration a M . To determine the Lagrange multiplier A, we solve the equation of motion for u M and obtain 
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2a"V a E{4>) - 4>E 



A< r 



(15) 



where A^ F ) = 2d\u ,1 V v F ufl . Note that, as expected, A = if the dispersion relation is linear, which is obvious since in 
that case the action (|J) does not depend on the vector field u^ L . 

To conclude this part, we derive the energy density and pressure as defined by an observer comoving with u^, 

p=u»u»T^, p=^±^T^, (16) 

which gives 

P = 2 + \ V a <^ a 4> + ]T b np V 2n <t>V 2 ^ + uWT$ + 2A< F > , (17) 

n,p 

3p = <j) 2 - ^\7 a (f>\7 a <j) - 3 b np T> 2n (j)V 2p (f> + (3V 2 <j) + a a V</>) E{<j>) + S7 a E\7 a cf> + E0- Tj$. (18) 



Although the previous expressions hold for a non-interacting scalar field, it is trivial to include a potential term in 
the above equations. We now turn to the particular case of the FLRW metric. 
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C. Friedmann-Lemaitre-Robertson- Walker spacetime 



From now on, we restrict ourselves to a FLRW universe with flat spatial sections, the metric of which is given by 
g^Ax^Ax 1 ' = — At 2 + <z 2 (t)£y -Ax 1 Ax 3 , where a(t) is the scale factor and t denotes the cosmic time. The Kronecker 
symbol Sij represents the metric of constant time hypersurfaces in cartesian coordinates. In this particular case the 
expansion of the universe provides a natural definition of space and time. We thus choose the scalar function q to be 
the cosmic time t so that u M = — 5°, it' 1 = 6q. As a consequence, we now have -F M „ = 0. Also, the projector _L Mt/ takes 
the simple form _L MI/ Ax^Ax" — a 2 (t)8ijAx z Ax J . Using the previous equations, it can be checked that 



(19) 



i.e. Tif/D^Sjs the four dimensional expression of the three dimensional Laplacian, as expected from the general 
argument (|A5|). Let us note that the previous arguments can be easily extended to a FLRW universe with non-flat 
spatial sections. In a FLRW spacetime the commutation relations mentioned above are trivially satisfied and the 
energy density and pressure reduce to 



1 v 



1 •, 



I l.p 



di<l>dj(t> + ^b np 



I l.p 



g(n + p)-l 



v 2n <t>v 2p <t>, 



(20) 
(21) 



where we have implicitly integrated by parts and discarded a total derivative term for the pressure. 

We now derive the field equation for <f> to exhibit the modified dispersion relation, and in particular to establish the 
link between the series contained in the Lagrangian £ cor and the Taylor expansion that defines the dispersion relation. 
From now on, we use conformal time r\ defined by dt = a(t)Ar] and the reduced field fi = a<f>. Varying the action with 
respect to \i gives 



(22) 



where a prime denotes a derivative with respect to rj. The last term of this equation comes from the Lagrangian C coi 
which encodes all information on trans-Planckian physics. Since this term is composed only of Laplacians, it gives rise 
to a series in momentum, once we shift to Fourier space. For plane wave solutions of the form x) = [i,k(r\)e lk x 
one obtains the equation 



(J-k = 0, 



with the modified dispersion relation 

cj 2 (k) =fc 2 + 2a 2 ^(-l)("+rf& 



np 



n , p 



1 2(n+p) 



a{rfi 



(23) 



(24) 



This dispersion relation can be rewritten in terms of physical pulsation tt; phys = cu/a and frequency k phys = k/a, as 

<Jk) = fc p 2 hys + 2^(-l)(^)6 np fc 2 ^). (25) 

n,p 

Thus our procedure allows us to study the stress-energy tensor of a quantum field with dispersion relation such that 
uj 2 is an analytic function of k 2 

phys " plrys 

We now go one step further and second-quantize the system. We indeed study the behavior of gravitational waves 
and density perturbations in the trans-Planckian regime, which can be reduced in the pertubative approximation to 
the study of a free minimally coupled quantum scalar field propagating in a FLRW classical background, see e.g. [Bj. 
We thus write the field operator /} (a hat denoting an operator) in the Heisenberg representation as 



Ak 



(2tt) 3 / 2 



/ti(»?)e* x c fe +^()j)e 



ik-x - 1 



(26) 
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where C& and cj„ represent annihilation and creation operators of quanta with momentum k respectively, normalized 
as usual through = S(k — p). The vacuum state |0) is defined accordingly by Cfc|0) =0 and a star denotes 

complex conjugation. It is then easy to obtain the energy density and pressure of the scalar field 4> in its vacuum 
state by inserting Eq. (j|^) into Eqs. (|o|) and (|2~l| ) to get 



(0\p\0) = 
(0\p\0) = 



1 



4ir 2 a 4 
1 

47r 2 a 4 



dkk 2 



dkk 2 



l (k)\fi k \ 



2 2 duj 2 

3 dP 



^ 2 ) l^fel 2 



(27) 
(28) 



The above is one of the main results of the present article, and it constitutes the basis of the discussion to follow in 
Section III . For the standard dispersion relation of massless quanta uj(k) = k, one recovers known results, namely that 
in the limit pk oc e~ lujri , corresponding to subhorizon modes (see Section III), p — p/3 corresponding to a relativistic 
fluid, while in the limit (pk/a)' ~ 0, corresponding to super-horizon modes (see Section III), p = —p/3. It is of 
interest for what follows to note that even for the standard dispersion relation, the fact that the quantity p is frozen 
on superhorizon scales does not imply that the corresponding equation of state is of the cosmological constant type 
p = -p. 

For a modified dispersion relation, as mentioned before, the energy density is the straightforward generalization 
of the standard expression, with uj 2 {k) denoting the modifed pulsation. However the pressure expression is a non- 
trivial generalization, and the presence of the term duS 2 /dk 2 implies that various equations of state may be obtained, 
depending on the shape of the dispersion relation (this latter also determines the mode function p k ). 

In the following, we will also be interested in the power spectrum of the fluctuations of the field 4>, defined by 



(o|^,x)|o) = y T k 3 P(k), 

where the power spectrum per logarithmic interval k 3 P(k) reads 



(29) 



1.3 

em - ±, 



Pk_ 

a 



(30) 



In the case of power law inflation and a standard dispersion relation, the power spectrum is given by k 3 P(k) 
Ask ns . In the particular case of a de Sitter spacetime, the spectrum is scale invariant, i.e. ng = 1. 



III. APPLICATIONS: EQUATION OF STATE AND BACKREACTION 



In this section we apply our previous calculations of the s tress-energy tensor to two different situations that have 
arisen recently in the literature. As a first application (§ [II A) we calculate the equation of state of gravitational waves 
of super-Planck frequencies, which have been proposed to account for the observed accele ration of th e Univ erse [po[ ; 
our conclusions are different and do not support the claim made by these authors. In § [II B and § [II O] we study 
the backreaction problem of trans-Planckian modes in inflationary cosmology Jl5|,|l6|]. Our calculation is indeed well 
suited to this case as it allows us to explicitly calculate the energy density of the fluctuations and compare it to the 
background energy density. 



A. Trans-Planckian and dark energy 



It was proposed recently that gravitational waves of super-Planck frequencies with a dispersion relation that ex- 
ponentially decreases beyond the cut-off (Planck) momentum would contribute significantly to the present energy 
density, with an equation of state mimicking a cosmological constant ppfl . To discuss this claim, let us consider the 
following dispersion relation 

= ^ = k 2 e- k ^\ (31) 

where k c is the cut-off momentum. Even though it is not exactly the dispersion relation considered in Ref. p0| , it 
is simpler and very similar: it is linear at small momenta (k <C k c ), reaches a maximum around k c and decreases 



G 



exponentially at high momenta (k 3> k c ), and it will serve our purposes well enough (in addition, it was claimed that 
the result does not depend on the details of the dispersion relation). 




FIG. 1. Dispersion relation with an exponential decreasing tail, as envisaged in Ref. |20|. The dot-dashed line shows the 
linear dispersion relation, and the thick solid curve shows the modified dispersion relation. The dotted lines delimit three regions 
of different evolutions for the mode function: in region I, corresponding to the "tail", a mode has trans-Planckian momentum 
k > ak c and sub-Hubble comoving pulsation u(k) < aH; in region II, a mode has ui(k) > aH , and in region III, modes are 
super-horizon sized k < aH. The dashed area represents the "tail" modes defined by k > K + (equivalently fc phys > K + /a in 
the figure). 

In Fig. [I] we show the qualitative behavior of the comoving pulsation associated to this dispersion relation as a 
function of comoving momentum and conformal time. There exists a region ("tail" modes) in which the modes have 
^phys > k c and a pulsation smaller than the Hubble frequency kj p hy S < H. These modes are defined by k > K+ > ak c , 
where K + is obtained by equating u> with the comoving Hubble rate aH, 



tf + ~afc c W21n^|. (32) 

We further assume that the scale factor evolves as a power law of conformal time, i.e. a{rj) = aoiv/vo)^ , where ao is 
the dimensionless scale factor at time 770, and (3 — 1,2,-1 respectively for radiation domination, matter domination 
universe and de Sitter inflation. It is important to note that a mode that is contained in the tail at time r\ was 
already contained in the tail at any previous time; indeed it can be easily verified that the wavenumber K + increases 
with time for an expanding scale factor and constant or decreasing curvature. In the "tail" region modes have a 
pulsation smaller than the Hubble expansion rate, i.e. u> <C aH — (3/r/, which implies ui 2 <C a" /a = f3(/3 — l)/^ 2 - 
As a consequence, the term uj 2 in the equation of motion (|23|) can be neglected with respect to a" /a, so that the 



approximate solution to Eq. (23) reads 



n dn' 

fx k ( V ) ~ C + (k)a( V ) + C_(fc)a(r?) J (33) 

where C+ and C_ are the coefficients of the growing and decaying modes respectively. These two functions are inter- 
related through the Wronskian normalization condition nul^X ~~ Vkt^'k ~ ^ ne solution is indeed "frozen" since for 
the growing mode the scalar field <p = /ife/a is constant in time. The claim of Ref. |2(J is based on the assumption that 
frozen modes have an equation of state of a cosmological constant type p = —p. As already mentioned above, this is 
not true even if the dispersion relation is standard and the mode is frozen (super- horizon sized), since the equation 
of state then reads p = — p/3. In addition, as we argue below, this claim does not hold either when the dispersion 
relation is modified. 
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In order to calculate the energy density and pressure, and the relation between these two quantities, it is necessary 
to specify the vacuum state. Unfortunately, for the "tail" modes, it cannot be chosen unambiguously since the WKB 
approximation breaks down in the limit cu <C aH and the concept of adiabatic vacuum cannot be used. This ambiguity 
also implies that the initial value of fj,, i.e. the coefficients C—(k) and C+(fc) cannot be determined unambiguously. 
This problem has not been addressed in Ref. as only one branch of the solution to the equation of motion is 
considered and the normalization coefficient (similar to the above C±) is taken to be independent of k. In any case the 
solution given in Ref. |2(J is not a correct solution to the field equation, since these authors neglect the contribution 
of the term a" /a in the limit r\ — > — oo, but a direct comparison between the various terms in their Eq. (22) shows 
that instead the term a" /a dominates in this limit, in agreement with the above discussion. 

In the following, we propose to circumvent the ambiguity related to the definition of a correct vacuum state by 
discussing two possible generic sets of initial conditions. 



1. Power-law initial conditions 



The above problem is in a certain sense similar to the situation encountered in cosmology before the advent of 
inflation. Indeed in the absence of an inflationary epoch, quantum modes can only enter the horizon as time goes, and 
the initial data had to be specified on scales larger than the horizon in a regime where the field is frozen. In contrast 
inflation has the virtue of stretching modes beyond the horizon, so that the initial data can be specified without 
ambiguity while the mode lies well inside the horizon. In the absence of such a mechanism, it had been proposed 
to use a universal power law to describe the power spectrum of the initial fluctuations, an approach that relied on 
astrophysical arguments. We thus propose to adopt a similar choice and take the growing mode coefficient to be a 
power law in momentum 

-1/2 ( k 



C + (A)=*r^l-j , (34) 

where kj is a constant related to the amplitude of the fluctuations and a a constant. It follows from Eqs. (|27|-|2^) that 
the energy density and pressure of the modes in the tail at a time r\ are 

<0|p|0> = ^ / u^e^^du, (35) 



K+/kj 



(0|p|0) = --^W / u4+2 ° [! + 2X2u2 ] e~»' x 'du, (36) 



K+/ki 



where we have introduced the dimcnsionlcss parameter A which depends on the time at which the density and pressure 
are evaluated A = ki/(ak c ). This can be integrated to give 

(afc c ) 4 (ak c \ 2a+ \( , 5 K\_\ 

!„/ 5 Ki \ 2 / 7 K 



(o|p|o>- {akc) (ak 



1 



<°I» = W f T [ a+ r^hh (37) 



87r 2 a 2 V h 



r u + - — t- +_r 



3 V 2' a 2 B 3 V 2' a?k? 



+ 



(38) 



where T(f3, x) = f°°dte H@ 1 is the incomplete gamma function |p4j| . The equation of state parameter is thus 
obtained as (note that it does not depend on the arbitrary wave number fci) 

= P 1 2 T [a + 7/2,^/^)] 

3 3T [ a + 5/2,K 2 + /(a 2 k 2 )] ' 1 ' 

and since K + ak c , one can use the large argument expansion of the incomplete gamma function T(P,x) ~ 
x^~ x exp (—a;) at large x to derive the following expression 

1 2 ;«-VftV (40) 



3 3 3 \H _ 

where we have used Eq. ( |32] ) to express K+/(ak c ). If k c is of order of the Planck scale, with Ho/Mpi ~ 10 -61 , then 
one obtains today 
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186. 



(41) 



This value does not depend on the normalization scale k\ and is in clear contradiction with a cosmological constant 
type equation of state, i.e. p ^= —p. It is rather a "phantom energy component" according to the terminology 
introduced in Rcf. po| , i.e. a component for which u> < — 1. In Rcf. p5| , observational constraints on the equation of 
state have been studied, and ui ~ —186 already appears in contradiction with the SNIa data, see Fig. 6 of Ref. [p5| . 
Furthermore, if one evaluates the energy density contained in the "tail" modes today (using again the asymptotic 
value of the gamma function), one finds 



Ptail 



H 2 k 2 



2a + l 



(42) 



or if we rather calculate retail = Ptaii/Pcrit, where p cr it = 3H 2 Mp l /8ir is the critical energy density, one obtains 



tail 



o.i 



Mi 



pi 



(43) 



Therefore, if k c ~ Mpi, f2taii ~ 1 if a = —1/2 and/or k\ — K + ~ 16.76Afpi. But the important point is that this 
can be obtained only at the expense of fine-tuning of the parameters, and the above choice docs not seem natural; 
actually, the parameter k\ is arbitrary and a priori there exists no natural way to fix it. 



2. Minimising energy 



A second possibility is to fix the initial conditions by requiring that the initial state configuration minimizes the 
energy, as advocated in Ref. |2(| and used in Refs. To this end we parametrize fj,' k as fJ-' k /pk = x + iy, where x 

and y are two real functions, and we rewrite the energy density as 
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(44) 



with H. = aH . One can now find the extrema of this expression under the Wronskian normalization constraint to 
obtain the mode function and its derivative at the initial time rji 
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It follows that in the "tail" region, the field /i evolves in time as 



VkiVi) = [W + iw(fc,f7i)]/i fc (77i). 



^2w(fc,77i) a(r)i) 



1 + lLu(k, T]i) 



0(77) 



dij 



The energy of the "tail" modes at a later time r\ can then be estimated as 
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(47) 



It can then be checked that this quantity depends strongly on the choice of the initial time. Moreover, since the 
dispersion relation decreases exponentially fast as k — > +00 and rji < 77, the integral diverges exponentially when 
k — > 00. 

In spite of the ambiguity related to the choice of the initial state, the above two solutions have the merit to show 
that there is no rigorous argument in favor of the claim made by Mersini et al. |po| . In particular, one does not 
generically find (or expect) the "frozen" modes of the "tail" neither to have a cosmological constant type equation of 
state nor to have an energy density today coinciding naturally with the critical energy density. As we have mentioned 
above, the difference with the results obtained in Ref. [pof presumably lies in an incorrect assumption made by these 
authors to derive the time evolution of the mode functions. One should also add that these authors |2(]] have not 
addressed the problem of the energy contained in modes with momentum k < K + (in particular region II in Fig. [l]). 
Indeed if one considers the "tail" modes as a source of gravitational energy today, there is a priori no reason to 
discard the contribution of other modes. However the subhorizon non-frozen modes oscillate, and it is easy to see 
that their energy density will be of order Mp h which is nothing less that the celebrated long-standing problem of the 
cosmological constant. 







B. Backreaction: general discussion 



28| (in collaboration 



The previous section has briefly touched upon the issue of backreaction in trans-Planckian cosmology, which arises 
whenever the energy density of the fluctuations (scalar or tensor) becomes comparable to the background energy 
density. In this case, the perturbative approximation used to describe the evolution of the mode functions breaks 
down, and one has to resort to a second order (complex) treatment of the Einstein equations p7 

Backreaction effects may arise in different situations in inflationary cosmology. Liddle and Lyth 
with E. Stewart) pointed out that if the modes of comoving wavelength of order of the horizon size, i.e. modes of 
cosmological interest for structure formation, start in a non-vacuum initial state, then inflation suffers from a very 
important backreaction problem. It is interesting for the following discussion to summarize briefly their argument. 
A mode with momentum k exits the horizon N c (k) ~ 56 — \n(k/Ho) e-folds before the end of inflation P9| . At 
horizon crossing, this mode carries a physical momentum fc p hy S = -ffinf , -ffinf ~ 10 13 GeV being the Hubble scale 
during inflation. Therefore, this mode had a physical momentum fc p h ys = exp(Ai)_ff; n f at the onset of inflation, 
N[ = iVtot — Ac being the number of e-folds of inflation before horizon crossing, and A to t the total number of e-folds 
of inflation. There is no reason to expect that inflation lasted just sufficiently long so as to stretch beyond the horizon 
only those modes with wavelength smaller than the horizon size today, and in general A tot ^S> A hence Ni 3> 1. The 
energy density contained in these modes at the begining of inflation is thus 5pk ~ (27r 2 ) _1 n/ c exp(4Ai)i?^ lf up to a 
numerical factor of order unity, with rifc the occupation number (we discard the zero-point energy). The ratio of this 
energy density to the critical energy density thus reads Sp/p cr i t ~ nu exp(4A;)(i?i n f/Mpi) 2 ; finally H in f ~ 10 -6 Mpi 
implies the very stringent constraint rik ^$ exp(— 4Ai + 28) to avoid backreaction problems. In other words, either 
./V; < 7, which may seem unnatural, or nt < 1 which implies that these modes are in their vacuum state or very close 
to it. Even though this argument was developed in the context of standard inflationary cosmology without modified 
dispersion relation and trans-Planckian physics, it can be carried over directly to our present problem. 

In fact, this has been done recently by Tanaka |15| who assumed that the modification of the dispersion relation 
could be represented by a non-vacuum initial state in a theory with unmodified dispersion relation. This implied, 
in much the same way as above, that backreaction effects should be very important. Starobinsky [ fl6[ has further 
argued that if adiabaticity is broken at some point in the evolution of a given mode, then at late times, when 
adiabaticity is restored, this mode behaves as the sum of positive and negative energy plane waves, meaning that a 
finite amount of quanta has been created, and the energy density they carry is large compared to the critical energy 
density today. However neither of these authors had at their disposal a rigorous expression for the energy density and 
they could not provide a quantitative estimate of this backreaction energy density. Nevertheless as we show in the 
following their conclusions remain correct, notably because the correct expression for the energy density Eq. ( p7j ) is 
the generalization of the standard expression. The overall picture for the origin and magnitude of backreaction can 
be described as follows. 

Consider a general dispersion relation such that for a given comoving wavenumber k, for time 772 (fc) < r/ < 773 (fc), 
where r]^{k) ~ —1/k denotes the conformal time of horizon crossing, the WKB approximation is valid, i.e. u>' /uj 2 <C 1. 
This notably implies cu(k) ^> H for i]2(k) < r\ < 773 (k), as the WKB condition would not be satisfied otherwise. In 
this case, the mode function evolves according to 



MfeO?) 



a(k) -i r w (fe,V)dV , /3(fc) i P u(k,v')dv' 



VMM) 



VMM) 



(48) 



with |ce(fc)| 2 — |/3(fc)| 2 = 1 from the Wronskian normalization condition. In principle the Bogoliubov coefficients a(k) 
and (3{k) depend on time, but their time dependence can be neglected to first order in the WKB expansion. These 
coefficients can be obtained by matching the mode function pk and its first derivative at time 772(A) with the solution 
of the field equation in the region 77 < 772 (fc) - If the WKB approximation was also valid at all times 7/ < 772 (fc), and the 
Bunch-Davies adiabatic vacuum is chosen as the initial state of the field (in mode k at least), then (3k ~ and a/- ~ 1 
at time 77. Furthermore, the modification to the power spectrum in this case is of first order in /3j., and thus remains 
small juj. However, if at some time prior to ^(fc), the WKB condition was violated, then |/3fc| can be a priori large; 
this had been remarked by Starobinsky [ flGf . This may induce large modifications in the power spectrum, but it also 
represents the creation of a large amount of energy density due to the breaking of adiabaticity. In effect, the energy 
density contained at time 77 in all modes such that the mode function (ik(v) is given by Eq. (|4£]) can be written as 
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1 



47r 2 a 4 



+ 



dfc/e 2 
2u> 



1 

2lu 



\0k\ 



|7| 5 



2uj 



2 1 2 
u + 7 



^ — 2i J 71 cj(fe ) ?7 / )d77 / 



(7* 



2i J ' a;(fe,7?')d?? / 



(49) 
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where we have used |a/c| 2 = 1 + \Pk\ 2 ■ Again, note that the above integral only covers the range of comoving momenta 
for which Eq. |l8| ) is valid. In principle, the total energy density should also contain the contribution of other domains 
of momenta. In the above expression, the quantity 7 is defined as follows 



7(M) 



Lj'(k,T]) ^\ _|_ a ' 



(50) 



In a situation where WKB is a good approximation, we have j/cj ~ i and the previous expression reduces to 



(0|p|0) 
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(51) 



Note that in order to remove the two oscillatory terms, no procedure of time averaging is needed in contrast with 
what has been done in Ref. If uj{k) does not decrease faster than fc -3 as k — ► +00, this integral diverges as in 
flat space due to the first term in the integrand. We interpret this infinite quantity as the zero-point energy which 
we subtract (adiabatic regularization), and interpret the remainder as the presence of finite energy density in modes 
with occupation number \(3k\ 2 , as usual. It is interesting to note that if u>(k) decreases faster than fc~ 3 as k — > +00, 
the integral is no longer divergent, and delicate questions on the necessity of renormalization arise. We will not touch 
upon these subtle issues in the present article, and will adopt the simplistic point of view in which the zero-point 
energy is subtracted in all cases. In Ref. Q, an expression similar to Eq. ( |5l] ) above had been used to discuss the 
magnitude of backreaction, but the origin of (3k had been left unspecified. The above discussion establishes the link 
that was missing and it also justifies more rigorously the approach of Refs. 15 R| , i.e. it shows explicitely that having 
a non standard dispersion relation is equivalent to considering non-vacuum quantum states for the perturbations as 
guessed in these studies. Furthermore, it also shows how to calculate (3k, i-e. through the matching with the solution 
to the mode equation at time 77 < r)2(k). 

In the following, we provide a concrete example of such a calculation and discuss the magnitude of (3k for a general 
class of dispersion relation. The WKB approximation is valid whenever |a//a; 2 | < 1, or 
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< 1. 



(52) 



The WKB approximation can thus be violated in two general ways: by space-time curvature effects, i.e. when 
Wphys ^ > and/or by singularities in the dispersion relation or its first derivative. In the following we will be 
interested only in the former class. An illustrative example of the latter is given by the Corley-Jacobson dispersion 
relation that becomes compl ex, w hich has already been discussed extensively We note that the dispersion 

relation considered in Section III A violates the WKB approximation at large physical momenta in both ways, namely 
d In w p hys/d In fc p hys diverges as fc p h ys — * +00 (or 77 — > —00), and w p h ys decreases much faster than H as 77 — > —00. In 
order to compute the Bogoliubov coefficient (3k one needs to solve the mode equation in the region where the WKB 
approximation is violated. However, as discussed in Section III A, in this region one cannot determine unambiguously 
a vacuum state, and in particular one cannot fix unambiguously initial data. We thus make a further assumption, 
and assume that the dispersion relation is such that the WKB approximation is restored in the far past 77 — * — 00. 
More precisely we assume that for a given comoving wavenumber fc, there exists a time 771 (fc) such that for 77 < 771 (fc), 
the WKB condition is satisfied; this allows us to define proper initial data for the mode evolution. 



C. Backreaction: examples 



A dispersion relation with the above property can be constructed easily. For instance, consider the lowest order 
generalization of the Corley-Jacobson dispersion relation, i.e. w 2 hy = fc 2 hys + 26nfc 4 hya — 2&i2fc® h s . Depending on the 
sign of the coefficients bn and 6 12 , this dispersion relation can present both a maximum around fc c and a minimum 
at physical momentum larger than fc c . For a specific choice of these coefficients, this minimum can be chosen to be 
smaller than the Hubble scale at some time 77, as depicted in Fig. ||[ 
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This is only an example of the desired class of dispersion relations, and the following discussion remains general 
with respect to the form of w p hys(fc p hys)- We fix the initial conditions in the region r\ < rji(k) where the WKB 
approximation holds by selecting positive frequency modes, which corresponds to the choice of the Bunch-Davies 
adiabatic vacuum state for the field, 
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At time r?i(fc), this mode enters region I where u>{k) <C H, and and its first derivative must be matched to the 
solution Eq. ([33|), which gives 
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dry' 



(54) 



where the quantity 7 has been defined previously, see Eq. (J50|). Note that 71 carries dimensions of inverse time and 
in particular, for general C 1 dispersion relations, 71771 is a number of order unity, which depends weakly on k [see 
Eq. (j(3(]) below]. At time T] 2 (k), the mode enters region II in which the WKB approximation holds, i.e. uj(k) Ti, 
and /ifc and its first derivative must be matched to the WKB solution, giving 
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and the two coefficients a k and 0k are explicitely given by 
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(56) 



(57) 



and where 72 is defined as in Eq. ( |50| ) albeit it is evaluated at time 772 (k). The previous equations give the general 
expressions of the coefficients a k and (3 k in region II. These expressions can be further worked out in the following 
two situations: (i) when |tji| ^S> 1 772 j and the scale factor is written as 0.(77) = a inf (v/Vinf Y or (ii) when 772 ~ 171. 
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Let us start with the first situation. Assuming \rji\ 3> |??2| and using W2772 = <^>iVi 
contribution to the two coefficients and (3k can be written as 
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The values of 71771 and 72% can be cast in the form 
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(60) 



which shows that these quantities are of order unity and do not depend strongly on k unless the dispersion relation 
contains exponential factors important at the time of matching. It follows that in region II, one has 



4%) 



1 



2/3-1 
V2\ 1 ,2 

72% 

771 



7i m 



1 - 2/3 



cos 



2^ + 2 / oj(k,r]')dri' 
m(k) 



(61) 



where </?2 = arg(72). In region III, in which the mode exits the horizon on the linear part of the dispersion relation, the 
solution reads fJ>k(v) — Cm(k)a(r]). The constant Cm{k) is obtained by matching the super-horizon solution with the 
solution in region II given above and it essentially determines the power spectrum [see Eq. (fjO])]. When the mode k 
exits the zone II and enters the zone III, the dispersion relation is linear w ~ k. We deduce that 773 (fe) = — y73(/3 — l)/k 
and thus that 
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Two interesting features should be noted on this form of the power spectrum. First, the standard spectral index of the 
overall amplitude ng — 1 = 2/3 + 2 is modified due to the factor (?72/?7i) 2 ^ 1 which a priori depends on k. Second, the 
power spectrum exhibits superimposed oscillations and therefore it can a priori vanish for some values of k. The above 
formula can be used to study the features of these oscillations once the dispersion relation has been specified but we 
will not pursue this goal here. However, in the particular case of a de Sitter inflationary period (/3 = — 1), the previous 
features are no longer present since all quantities are functions of kr\ only. Notably, since H is a constant in de Sitter 
space, the two solutions 771 and 772 to the equation w p h ys = H read 771 — C\jk and 772 = C^/fc, where C\ and C2 are 
constants, and therefore 772/771 is independent of k. For the same reason, since u>(k,rj) = u)(krj), a change of variables 
i] — > u = kf] in the integral appearing in Eq. ( |62| ) shows that this quantity does not depend on k. The spectral index 
is thus unchanged, i.e. ng = 1, and the superimposed oscillations reduce to a constant numerical factor. The power 
spectrum is thus unchanged except for an overall modification of its amplitude; this latter is actually magnified by a 
factor O [(771/7/2) 3 ] (/3 = —1 for de Sitter). This result is in agreement with Refs. (l^,[l3|] where the power spectrum 
was calculated for a theory with modified canonical commutation relations but unmodified dispersion relation. In 
Ref. 12[ , it was argued that the power spectrum spectral index remains unchanged for de Sitter inflation, and in 
Rcf. [ 13| it was further argued that only the overall amplitude of the power spectrum is affected. 
Finally, the energy density contained in the modes in region II at time 77 can be written as: 
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and the domain of integration K, is such that aH < k < aki with fc 2 being the smallest wavenumber such that 
Wphysfe) = H and k 2 > H . It is artificial to consider different domains of wavenumbers and to treat the corresponding 
contributions to the energy density separately, but those modes with wavenumbers k > ak\ [where k± is the other 
wavenumber > H such that cj p hy S (fci) = H] are in their vacuum state and WKB holds so the contribution to the 
integral has been removed, while those with ak\ < k < aki are in the non-WKB zone, where it is ambiguous to define 
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a vacuum state and to calculate the contribution to the energy density. On the other hand, it is sufficient to show that 
one of these contributions is of the order of the background energy density to demonstrate that there is a backreaction 
problem. This is the spirit of the following calculation and the interval chosen is particulary well-suited since in the 
standard case the corresponding contribution vanishes. In the present case, one obtains for the case ris = 1 



/ dk phys kl hys uj phys (k) = ( — ) 0(k*), 
Jk V72/ 



(64) 



where we have approximated the integral to the peak value of w p hys at k c . This value should be compared with the 
background energy density during inflation pi n f = Al^H 2 ^. Therefore, if we take k c ~ Af P i, we see that p 3> Pi n f due 
to f/i n f < Mpi and 771/772 1, the latter equation expressing the fact that the coefficient (3k is large in region II. This 
result coincides with the analysis of Refs. |l5| , |l6| ]. 

The other situation where the spectrum and the energy density can be evaluated explicitely if the time 772 is close 
to 771. If we write 772 =771(1 + e), where e is a small parameter, then one finds for the two coefficients a% and (3k [see 
Eqs. <M,E%] 
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where po and po are the background energy density and pressure at time 771 (fc) respectively. The coefficient Q is 
defined by Q = 3o/ 2 /(4w 2 ) — lj"/(2lj). The interpretation of these formulae is as follows. To leading order in e, the 
coefficient oik is just a pure phase and the coefficient (3k vanishes as expected. Two terms appear to the next order in e 
where the corrections show up. The first is Q/oj 2 . The condition for WKB to be valid is \Q/oj 2 \ <C 1 and therefore this 
term indicates the magnitude of violation of the WKB approximation. The second term is proportional to po/3 — po- 
The presence of this term is also natural because it vanishes for radiation (a" = 0), in which case the exact solution 
for the mode function Pk(jl) is a complex exponential and no correction is expected for the overall amplitude. Finally, 
repeating the same calculations as above for the power spectrum, one finds 
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(67) 



The correction to the power spectrum is of order e as expected. To this order, as before, we have a modified overall 
amplitude and superimposed oscillations appear. We can also estimate the energy density. Inserting the expression 
of the coefficient (3k into Eq. (|5l|), one obtains 



P - 0(e 2 k 4 c ). 



(68) 



Therefore, there is no backreaction problem if e 2 k 4 < M^, x Hf n{ . Let us write k c as k c = 10~ s Afpi, where the coefficient 
s fixes the scale of the characteristic wavenumber with respect to the Planck mass. One can also write e as e = 10 -p , 
where p roughly gives the order of magnitude of the modification to the power spectrum. If we take 7J inf = KT 6 Af P1 , 
then there is no backreaction problem if p + 2s > 6 that is to say if s > 3 — p/2. Therefore a modification of the 
spectrum of order 10% in principle already detectable now with COBE or in the near future with MAP can exist 
without a significant backreaction problem provided the characteristic scale k c < 10 _2 5 Mpi (but larger than the 
Hubble scale of inflation). Similarly, a modification of 1% in principle observable by the Planck satellite mission can 
be obtained if s > 2. Interestingly, these domains encompass the Grand Unification scale ~ 10 16 GeV and possibly 
the string scale. 

At this stage one should note that the issue of backreaction in the present Universe []l6| , which can be interpreted 
as the production of gravitons of super-Planck momentum, applies to those dispersion relations which violate WKB 
today, i.e. those for which u> < ciqHo today (notwithstanding singularities in the dispersion relation). Since the 
comoving Hubble scale today is orders of magnitude below the comoving Hubble scale of inflation, a dispersion 
relation which broke the WKB approximation during inflation does not necessarily imply production of gravitons 
today. This holds in particular for those dispersion relations above which entail a modification of the power spectrum 
without a strong backreaction problem at the time of inflation. 
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Finally one should note that the above situation may be encountered for a wider class of dispersion relations, if 
the early time evolution of the scale factor is different from inflationary expansion. Notably consider a spacetime 
which is asymptotically Minkowski as 77 — > —00, with a scale factor evolving as a{rf) = a\ + a m i (77/77;^ )^ ■ Provided 
dlnwphys/dlnfcphys is not singular, the WKB approximation is valid asymptotically. One could reproduce the above 
calculation with the substitution for the new scale factor, and the final expression for the power spectrum would be 
the same [at late times, Oi becomes negligible compared to the expansion term in a{rf)\. However the ratio 772/771 then 
depends on k if 772 is in the far past when the term a; cannot be neglected even if [3 = — 1. 

IV. CONCLUSIONS 

We have derived the stress-energy tensor for a free scalar field with a general modified dispersion relation. In 
particular, we have obtained the expression of the energy density and pressure in a FLRW background. This result 
generalises previous studies jl|,[l{| which were restricted to the Corley-Jacobson dispersion relation. We have applied 
our calculation of the stress-energy tensor to a series of examples to discuss the equation of state of trans-Planckian 
modes and the issue of backreaction in inflationary cosmology. 

We have first examined in details the possibility that the energy contained in trans-Planckian modes with a frequency 
much smaller than the Hubble expansion rate could account for the form of vacuum energy density measured in the 
low-redshift Universe |2(i(| . In the case of dispersion relations with ultralow frequencies at high momenta, as proposed 
in Ref. p0[ , one cannot select a vacuum initial state without ambiguity. Nevertheless, using well- motivated proposals 
for this initial state, we have shown that the equation of state of these trans-Planckian modes does not have the 
correct form, contrary to the claim made in Ref. Moreover, the numerical value of the energy density is not of 
the order of the critical energy density unless fine-tuning is required, and we thus conclude that the scenario proposed 
in Ref. |2(J does not stand up to scrutiny. 

We have also discussed the issue of backreaction in trans-Planckian inflationary cosmology. In particular, we have 
focused on a class of dispersion relations for which a; p h V s < H for a finite time interval for trans-Planckian comoving 
momenta. This is the general class of dispersion relations for which the WKB approximation is broken at some point 
during inflation, but is restored in the far past as 77 — > —00. If the evolution is adiabatic all throughout inflation up 
to horizon exit, it is known that the power spectrum is not modified. In the above case, the WKB approximation 
is precisely broken in the region where w p h ys < H, but is valid at earlier and at later times. We have obtained the 
analytical expression for the amount of energy density stored in modes at late time, and find that it is in general 
much larger than the background energy density. We have computed the power spectrum of metric fluctuations and 
showed that this power spectrum is not modified except for its overall amplitude, independently of the dispersion 
relation and whether WKB holds or not, if the inflationary period is strictly de Sitter. This supports the belief that 
inflation is robust to a change in the dispersion relation. If however the inflationary period is not strictly de Sitter, 
then the power spectrum is tilted with respect to the standard case (unmodified dispersion relation) and superimposed 
oscillations appear. Finally we have exhibited a class of dispersion relations for which the power spectrum of metric 
fluctuations is strongly modified and the initial conditions can be set up properly, since they are fixed in a region 
where the WKB approximation holds. 

Our work thus completes the arguments developed in previous works |^,|3|,f^,[l4| on the relation between adiabaticity 
of the mode evolution, the modification of the predictions of inflation, and the magnitude of backreaction. In particular 
the following picture seems to emerge: if the evolution of the modes is adiabatic all throughout inflation up to horizon 
exit, then the power spectrum is unmodified (or weakly modified), and backreaction is not an issue. If however 
adiabaticity is broken at some point, then modifications to the power spectrum are likely to appear (except if the 
background spacetime is very close to de Sitter), but effects of backreaction then appear generic, and one must use 
a higher order expansion of the Einstein equations to derive meaningful conclusions. Finally, there exists dispersion 
relations such that the backreaction is weak but modifications to the power spectrum are not negligible, notably 
when the ratio k c /H- m f is not too large, when the time interval in which adiabaticity is broken is small, and when 
the background spacetime is not de Sitter. However this obviously requires fine-tuning, and overall a scale invariant 
power spectrum does indeed appear robust against changes in the dispersion relation if backreaction can be neglected. 
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APPENDIX A: NOTATIONS 



This appendix provides some identities satisfied by the vector field u^, the projector _L Mt/ and the derivative T> a that 
are implicitly used in the calculations of Section |ll|. Since the norm of is conserved, one has trivially UaV^u" = 0. 
The covariant derivative of m m can be conveniently decomposed as 

V M u„ = ^9 _L M „ +er( M „) + w^„] - a^u v , (Al) 

where A^^ = {A^ v + A v ^)/2 and A^^ = (A^ — A u ^)/2 are the symmetrised and antisymmetrised parts of A^ 
respectively. The shear a^ v and vorticity uj[^ v \ are tracefree and purely spatial, i.e. 

<W" = 0, cr^=0, u^g^ = 0, uvu" = 0, (A2) 

and et M = m^V^u^ is the acceleration of the observer comoving with u^. The quantity et M is spatial, i.e. a M u M = and 
one has a M = for a geodesic. Using these relations, the expansion rate 9 can be written as 



1G 



(A3) 



The integrability condition of the hypersurface £ implies that ui^ = which is indeed satisfied if the definition of it p 
given in the main text is imposed. The projection operator _L M „ verifies 



T£= 3, ^ u" = 0, 



I P I 7 — I 7 



(A4) 



For the FLRW case, it follows that = 3H, a M = cfyi, = cj m „ = 0. 

Finally, if we denote by a 1 the internal coordinates on the hypersurface £ defined by the embedding x p = ^(cr 1 ), 
then the three dimensional spatial metric is given by -Lij— g i _ LV diX^djX v and it follows that T> 2 (j) is given directly by 



v-L 



(A5) 



where _L is the determinant of Hence V^V^cf) is the three dimensional Laplacian as defined by the observer 



comoving with vP. With the decomposition Eq. ( Al), it is easily shown using the identity u a u /3 V a V pep = <p — tt a V a ( 
that 



V 2 6 = Uc 



a a \7 Q <, 



(A6) 



where □ = V^V* denotes the four dimensional d'Alcmbertian. 



APPENDIX B: DERIVATION OF THE STRESS-ENERGY TENSOR IN THE PARTICULAR CASE OF 

THE CORLEY-JACOBSON DISPERSION RELATION 



We detail in this appendix the simple case where only b\\ does not vanish in the corrective Lagrangian. This 
Lagrangian contains derivatives of the metric that must be varied to obtain the stress-energy tensor. With 
5^/ zr g/8g pl/ = —\/— 55^/2, this stress-energy tensor can be written as 



SC 2 



sc 



so, 



p9iJ.ii 



For the present case £ cor = — &n(2? 2 </>) 2 , it leads to 

1 



T^ v = d»<$>d v 4> - - (d a (t)d a (j)) g pv + 2\u p u v - b n (2? 2 0) V 



+bn < Wp 



.5 V 2 



Sg 



.5 V 2 



5d p g^ 



(Bl) 



(B2) 



The first term in brackets can be written, using Eq 
5 (2?V 



while the second term gives 
5(V 2 



V^V^ + 0V (m u v) + 0u (i y v) <l> + 2u a u (i y V) V a (f>- i_ a/3 r^V^ - j>g a,3 T a(3 ^u v) , 



1 r 



Sdpg^ 2 
The derivative of this equation finally reads 



2 _lp(m y^)0_ yp^ + 2$gP^u^ - $g pv u p 



(B3) 



(B4) 



--do 



.5 V 2 



-g 



Sdpg^ 



(B5) 



It can be checked that the terms involving Christoffel symbols in Eqs. (B3) and (B4) strictly vanish, a necessary 
condition to the covariance of the stress-energy tensor. Finally, we end up with 
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-bn (V 2 



2U<j) ± lu , +2 (J> + 0$) 9tnJ - (V 2 <j>) 
+Au( ll V a u v - ) V a (j> - 4:U (jJ ,V v )U a 'V a (f) - 4o( /1 V„ 
-2b n 



w ( ^ v) (v 2 4>) + 2 (v 2 4>) ■ u {l y v) <t>- ±^ W a cf>W a (v 2 < 



The expression of the Lagrange multiplier A results from the equation of motion for 

— 6C 



ou^ \J—g 



-9 



which gives 



A = fen 



(V 2 <f) U + 0<p- 2a a V Q c/>) -(V 2 <j>)<j> + A (F) . 



(B6) 



(B7) 



(B8) 



It follows that the energy density and the pressure measured by the observer comoving with it' 1 (i.e. p = T^u^u 11 
and p = T^v \} LV /3) are respectively given by 



P = 4> 2 + ~V a 0V a + 6n (V 2 cpf + uWTW + 2A« 
P=l^ 2 - ^V Q 0V> + b f (V 2 ^ [3 (P 2 0) + 2a«V Q 0] + (p\ 
+ 2 -^Lv a ^{V 2 ^)- l -^T$. 



(B9) 



(BIO) 



In the case of the Minkowski metric, the mean values of the previous expressions calculated in a thermal state reduce 
to the formulas obtained in Refs. ||Tl||l9|| . 

Finally, the field equation for (f> is obtained by varying the action to obtain 



-a ' 



50 



hence 



□0 = -2fen [V 2 (V 2 (j)) + 2a w V u {V 2 <$) + (V 2 (f)] 



(Bll) 



(B12) 



It can be checked that the general formulae given Section ||, when applied to the special case of the present Corley- 
Jacobson dispersion relation, reduce to the equations given in this Appendix. 
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